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ABSTRACT 


An  experimental  study  was  conducted  to  determine  the  drag  coef¬ 
ficient  of  inert  spherical  particles  accelerating  in  a  laminar,  non-reacting, 
incompressible  continuum  flow.  The  Reynolds  number  range  which  was 
covered  in  the  study  was  from  150  to  1700,  and  particle  sizes  ranged  from 
150  ju  to  450  j El. 

The  convective  flow  behind  the  shock  wave  in  a  shock  tube  was  used 
to  accelerate  the  particles.  The  particle’s  diameter  and  the  displacement 
versus  time  measurements  were  obtained  using  a  rotating  drum  camera 
in  conjunction  with  an  oscillating  light  source.  The  photographic  data, 
the  particle  density,  the  shock  speed,  and  the  initial  pressure  and  tempera¬ 
ture  in  conjunction  with  the  normal  shock  relations  were  combined  to  cal¬ 
culate  the  drag  coefficient. 

The  drag  coefficient  is  usually  considered  to  be  a  function  only  of 
Reynolds  number  and  acceleration  modulus,  however,  CD  varies  consider¬ 
ably  because  of  particle  roughness.  Experiments  with  HP  295  ball  powder, 
whose  surface  is  relatively  rough,  produced  results  which  were  as  much 
as  85  per  cent  higher  than  the  steady  state  curve,  with  the  increase  depen¬ 
dent  upon  the  relative  Mach  number  of  the  flow  about  the  particle.  Similar 
drag  coefficient  experiments  with  smooth  sapphire  balls  did  not  produce 
the  scatter,  the  higher  values,  nor  the  dependence  on  relative  Mach  number. 


xiii 


I.  INTRODUCTION 


1. 1  ORIGIN  OF  THE  PROBLEM 

In  order  to  improve  the  density  and  specific  impulse  characteristics 
of  solid  propellant  rocket  fuels,  small  metal  particles  are  mixed  into  the 
propellant.  During  combustion,  condensed  liquid  and  solid  metal  oxide 
particles  are  formed  from  the  combustion  products.  These  particles  com¬ 
prise  30-40  per  cent  by  weight  of  the  combustion  products  in  current  solid 
propellant  rocket  motors.  Because  of  inertia  effects,  these  particles  leave 
the  nozzle  at  lower  velocity  which  means  a  loss  of  momentum  and  hence 
a  loss  in  specific  impulse.  In  order  to  calculate  the  loss  in  specific  im¬ 
pulse  due  to  the  velocity  and  temperature  lags,  the  drag  coefficient  of  the 
particles  is  one  of  the  variables  which  must  be  known. 

Presently  most  of  the  specific  impulse  loss  calculations  are  made 
using  the  ’’standard  drag  coefficient  curve”  for  spheres.  This  curve  is 
only  valid  for  a  single  smooth  sphere  in  a  steady,  incompressible,  laminar, 
non-reacting,  and  continuum  flow  field,  conditions  which  are  certainly  not 
satisfied  in  a  rocket  nozzle.  Thus,  in  a  rocket  engine  a  particle  is  moving 
in  an  accelerating,  turbulent,  compressible  stream  whose  temperature  is 
different  than  that  of  the  particle.  The  particle  also  moves  from  the  con¬ 
tinuum  regime  to  the  slip  flow  regime  and  possibly  to  the  free  molecule 
regime  depending  on  its  size  and  the  flow  in  the  rocket  nozzle.  Figure  1 
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Figure  1.  Mach  Number — Reynolds  Number  Flow  Regimes  Encountered  by  a  Five-Micron 

Particle  in  a  Rocket  Nozzle 


represents  the  path  in  the  Mach  number  Reynolds  number  regime  which  a 
five-micron  particle,  produced  under  chamber  conditions  typical  of  solid 
propellant  rocket  motors,  can  experience^  \  In  some  cases  enough  par¬ 
ticles  are  present  such  that  the  particles  and  the  gas  must  be  analyzed  as 
a  two-phase  flow.  The  particles  may  remain  hot  enough,  due  to  the  tem¬ 
perature  lag,  to  emit  electrons  by  means  of  thermionic  emission  or  to 
change  the  drag  due  to  heat  transfer.  Some  of  the  particles  which  come 
through  the  nozzle  are  the  original  metal  particles  which  are  mixed  into 
the  propellant.  The  drag  coefficient  of  these  particles  will  be  a  function 
of  all  the  above  effects  plus  burning. 

In  order  then  to  determine  the  correct  drag  force  on  a  particle  in  a 
solid  rocket  motor  nozzle,  the  influence  of  all  of  the  above  parameters 
must  be  studied. 

1.2  REVIEW  OF  STUDIES  ON  THE  DRAG  COEFFICIENT  OF  PARTICLES 

Some  theoretical  and  experimental  work  has  been  done  to  determine 
the  drag  coefficients  of  particles  accelerating  in  the  incompressible  con¬ 
tinuum  flow  regime.  The  results  of  these  studies  differ  appreciably  from 
one  another  as  may  be  seen  in  Fig.  2. 

Ingebo^  conducted  experimental  studies  on  the  vaporization  rates 
and  drag  coefficients  for  isotane  sprays  accelerating  in  turbulent  air 
streams.  He  injected  the  liquid  drops  into  air  streams  moving  at  140 
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Figure  2.  Data  on  the  Drag  Coefficient  of  Spheres 


and  180  ft/  sec.  A  specially  designed  camera  was  used  to  obtain  drop-size 
distributions  and  drop-velocity  data.  He  was  then  able  to  obtain  vaporiza¬ 
tion  rates  and  drag  coefficients  for  the  liquid  drops.  He  found  his  results 
could  be  represented  by  a  single  curve  given  by 

27 

C  =  — — — 

D  Re’84 

Fledderman  and  Hanson'  ;  performed  similar  experiments  by  photographing 
spray  droplets  accelerating  in  streams  moving  from  50  to  75  ft/  sec.  Their 
results  are  one  hundredth  the  steady  state  value  for  a  sphere. 

Experimental  drag  coefficient  studies  of  burning  kerosene  drops  have 
been  done  by  Bolt  and  Wolf^  for  Re  <  1.  These  results  indicate  a  de¬ 
crease  in  the  drag  coefficient  due  to  burning. 

Rabin  et  al.  ^  determined  the  drag  coefficient  of  burning  and  non¬ 
burning  liquid  fuel  droplets  accelerating  due  to  the  convective  flow  behind 
a  shock  wave.  For  Re  >  200  Rabin's  drag  coefficients  were  greater 
than  those  of  a  sphere  in  steady  flow.  Rabin's  data  also  indicates  a  de¬ 
crease  in  drag  coefficient  due  to  burning. 

Rudinger^  also  used  the  convective  flow  behind  a  shock  wave  to 
determine  the  drag  coefficient  of  accelerating  glass  beads,  which  had  an 
average  diameter  of  29  microns.  The  x  versus  t  motion  was  recorded 
by  streak  photography.  The  Reynolds  number  range  of  the  experiment 
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was  from  40  to  300.  Rudinger  found  that  all  his  data  could  be  correlated 
by  the  expression 

6000 

D  _  1.7 
Re 

This  relationship  deviates  widely  from  the  steady  state  curve  for  spheres 
and  Rudinger  suspects  that  electric  charges  on  the  particles  may  be  the 
cause  of  the  deviation. 

A  very  thorough  literature  survey  on  "The  Fundamental  Aspects  of 
Solids-Gas  Flow"  was  made  by  Torobin  and  Gauvin^  ^ \  They  listed 
and  discussed  such  problems  as  the  sphere  wake  in  laminar  fluids;  ac¬ 
celerated  motion  of  a  particle  in  a  fluid;  the  effects  of  particle  rotation; 
roughness  and  shape;  and  the  effect  of  fluid  turbulence  on  the  particle  drag 
coefficient.  They  also  made  an  experimental  study  of  the  drag  coefficients 
of  single  spheres  moving  in  steady  and  accelerated  motion  in  a  turbulent 
fluid.  Small  radioactive  smooth  spheres  were  fired  into  a  turbulent  flow 
wind  tunnel  of  known  turbulence  intensity  and  the  motion  of  the  sphere  was 
recorded  by  means  of  a  radioactive  sensing  device.  Using  this  technique, 
continuous  x  versus  t  data  could  be  obtained  for  the  spheres.  It  was  found 
that  by  increasing  the  turbulent  intensity  level,  the  critical  Reynolds 
number  could  be  shifted  to  Reynolds  numbers  as  low  as  400  and  that 
drag  coefficient  was  independent  of  acceleration  in  turbulent  flow. 
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Recently  Crowe did  both  an  analytical  and  experimental  study  to 
determine  the  effects  of  burning,  evaporation,  and  acceleration  on  the  drag 
coefficients  of  particles  accelerating  in  gas  streams.  The  Reynolds  num¬ 
ber  range  extended  from  250  to  1600.  Both  the  analytical  and  experimental 
portions  of  the  study  were  confined  to  the  incompressible  continuum  flow 
regime. 

For  the  analytical  study  Crowe  chose  a  spherical  model  with  mass 
flux  through  the  surface  to  simulate  burning  and  evaporation.  The  tangential 
equation  of  motion  was  used  as  the  governing  equation  in  the  analysis  of  the 
boundary  layer  flow  about  the  sphere.  The  velocity  distribution  outside 
the  boundary  layer  was  assumed  to  be  that  corresponding  to  inviscid  flow. 
With  the  proper  boundary  conditions, the  equations  were  solved  indicating 
that  for  burning  or  evaporating  particles  the  skin  friction  coefficient  is 
reduced.  A  constant  form  drag  coefficient  was  assumed  and  used  for  all 
the  results.  The  total  drag  coefficient  is  then  just  the  sum  of  the  skin 
friction  and  form  drag  coefficients. 

For  the  experimental  portion  of  his  study,  Crowe  used  the  convective 
flow  behind  a  shock  wave  to  accelerate  the  particles.  A  high  speed  framing 
camera  was  used  to  record  the  particle  diameters  and  x  versus  t  history. 
With  the  shock  speed,  the  particle  density,  and  the  local  temperature  and 
pressure  known,  Crowe  was  able  to  calculate  the  drag  coefficient.  Crowe’s 
particle  size  varied  from  100  to  350  microns. 
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Crowe's  best  fit  curve  to  his  experimental  non -burning  data  is  approxi¬ 
mately  15  per  cent  higher  than  the  steady  drag  coefficient  curve.  His  ana¬ 
lytical  results  are  another  15  per  cent  above  his  experimental  data.  How¬ 
ever,  each  experimental  point  has  approximately  a  40  per  cent  probable 
error. 

1.  3  PURPOSE  OF  THIS  STUDY 

The  present  study  is  part  of  an  overall  investigation  to  study  the 
dynamics  of  inert,  reacting,  and  charged  particles  in  solid  rocket  motor 
nozzles,  and  is  a  continuation  of  Crowe's  work.  Basically  it  is  desired  to 
experimentally  determine  the  relation  between  C^  and  the  prime  variables 
Re,  M^/  VRe  ,  M^,  and  a  non-dimensional  burning  rate  parameter.  In 
addition  the  effect  of  the  secondary  variables,  relative  turbulent  intensity, 
unsteady  effects  in  the  particles  wake,  particle  rotation,  and  particle 
roughness  are  also  to  be  studied. 

Due  to  the  relatively  large  experimental  error  and  scatter  in  Crowe's 
data  which  would  obscure  slip  and  compressibility  effects,  it  was  necessary 
to  repeat  the  incompressible  continuum  flow  regime  before  conducting 
experiments  in  the  slip  flow  regime  with  compressibility  effects.  In  order 
to  determine  C^  more  accurately,  a  new  and  better  instrumented  shock 
tube  was  constructed.  This  study  is  primarily  concerned  with  a  more 
accurate  determination  of  the  drag  coefficient  of  small  spherical  particles 
accelerating  in  the  incompressible  continuum  regime. 
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H.  PARTICLE  EQUATIONS  OF  MOTION 

Newton’s  second  law  of  motion  states  that  the  rate  of  change  of  mo¬ 
mentum  of  a  particle  is  equal  to  the  sum  of  the  forces  which  act  on  the 
particle  and  is  in  the  direction  in  which  the  sum  of  the  forces  acts.  In 
mathematical  form  Newton's  law  is 

^F  =  mfl  (1) 

The  forces  are  the  viscous  and  pressure  forces  which  act  on  the  particle 
surface  and  the  body  forces  which  act  on  the  particle  mass.  If  the  pres¬ 
sure  and  viscous  forces  are  expressed  in  terms  of  a  drag  coefficient, 
then  Eq.  (1)  can  be  written  as 

cDf  I  UR|  UR  A+ mf  =  ma  (2) 

where  CD  =  drag  coefficient  of  the  particle 

a  =  acceleration  of  the  particle  with  respect  to  inertial  space 
p  =  density  of  the  fluid 

UD  =  relative  velocity  between  the  particle  and  the  fluid 
R 

A  =  projected  characteristic  area  based  on  particle  diameter 
m  =  mass  of  particle 
f  =  body  force  per  unit  mass. 

For  the  case  of  burning  particles  Eq.  (2)  would  contain  an  additional  term 
for  the  momentum  flux  from  the  particle’s  surface.  The  body  force  term 
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in  this  analysis  is  the  gravitational  force  term  that  acts  on  the  body.  In 
the  present  study  the  particle  is  accelerated  by  the  convective  flow  behind 
a  shock  wave  so  that  the  acceleration  due  to  the  viscous  and  pressure  forces 
is  much  greater  than  that  due  to  gravity,  i.  e. , 

a 

—  >  100 
8  ~ 

Thus  the  gravitational  force  term  will  be  neglected.  Equation  (2)  may  be 
written  as 


CD2  '®r'  5RA=« 


(3) 


Since  the  flow  velocity  and  the  acceleration  vector  are  in  the  same  direc¬ 
tion,  the  vector  notation  may  be  dropped  and  Eq.  (3)  becomes 


2m  a 

puR2A 


For  spherical  particles  of  uniform  density, 


4 p  ad 
P 


3PUr2 


(4) 


(5) 


where  p  =  particle  density 

d  =  particle  diameter. 
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Dimensional  analysis  indicates  that  the  other  similarity  parameters 
which  are  needed  to  properly  determine  the  drag  coefficient  of  a  non¬ 
burning  smooth  spherical  particle  in  a  laminar  flow  field  subject  to  com¬ 
pressibility  and  non-continuum  effects  are  Reynolds  number,  Re,  Mach 
number,  M„,  Mn/VRe,  and  acceleration  modulus,  A  .  In  terms  of  the 
present  notation  Reynolds  number  and  Mach  number  are 


Re 

mr 


pu^d 


a 


(6) 

(7) 


The  particle  flow  regime  is  determined  by  the  value  of  MR/  VRe.  The 

(14) 

different  flow  regimes  are  defined  as  follows  : 


R  -2 

__fr  <  10  continuum  regime 

VRe  ~ 

-2  Mr  -1 

10  <  — —  <  10  slip  flow  regime 

VRe  - 


10 


-1 


< 


i/rT  - 


transition  regime 


“r 

Re 


>  3 


free  molecular  flow  regime 
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As  seen  in  Fig.  1,  a  five-micron  particle  in  a  typical  solid  propellent 
rocket  nozzle  experiences  all  of  the  above  flow  regimes  as  it  travels  through 
the  nozzle.  Equations  (5)  and  (6),  along  with  experimental  data,  can  be 
used  to  calculate  the  and  Re  of  a  non-reacting  spherical  particle.  Al¬ 
though  this  report  is  concerned  only  with  the  incompressible  continuum 
regime,  the  calculation  of  M_,,  and  M_/VRe  are  needed  to  insure  that 
the  particle  is  in  the  desired  flow  regime. 
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HI.  EXPERIMENTAL  APPARATUS 


3. 1  CHOICE  OF  EXPERIMENTAL  FACILITY 

In  order  to  study  the  dynamics  of  solid  particles  in  rocket  nozzles, 
an  experimental  facility  must  be  able  to  produce  a  particle  environment 
under  which  the  following  parameters  can  be  studied: 

1.  Acceleration  modulus 

2.  Mach  number 

3.  Mach  number/  VReynolds  number 

4.  Burning  rate  parameter 

5.  Electric  charges  on  particles 

A  shock  tube  was  chosen  because  by  using  the  convective  flow  behind 
the  shock  front  in  connection  with  the  small  spherical  particles  (50-500 
microns  in  diameter),  it  is  possible  to  produce  the  flow  conditions  under 
which  the  above  parameters  can  be  studied.  This  type  of  facility  was  used 

(13) 

successfully  by  both  Crowe v  '  in  studying  the  drag  coefficient  of  solid 

(5) 

particles  and  by  Rabin'  7  in  studying  the  shattering  of  liquid  drops. 

3.2  DESCRIPTION  OF  FACILITY 

The  experimental  equipment  consists  of  a  horizontal  shock  tube  into 
which  the  particles  are  injected  and  then  accelerated  by  the  convective 
flow  field  behind  the  shock  front.  The  distance,  x,  versus  time,  t,  of  the 
particles  is  recorded  by  a  modified  AVCO  rotating  drum  camera  in  conjunction 
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with  a  high  voltage  switching  circuit  which  supplies  energy  to  a  Xenon 
flash  tube.  By  combining  the  x  versus  t  data,  the  shock  strength,  and 
the  initial  conditions,  the  drag  coefficient  of  the  particle  can  be  obtained. 

3.  2-1  Shock  Tube  and  Driver  Section 

The  13/8  inch  square  shock  tube  consists  of  a  three  foot  stainless 
steel  driver,  a  six  foot  stainless  section  between  the  driver  and  the  test 
section,  a  one  and  one -half  foot  aluminum  test  section,  and  a  six  foot 
stainless  section  downstream  of  the  test  section. 

The  diaphragm  material  is  ruptured  by  means  of  a  long  rod  which  is 
inside  the  driver  section  and  is  driven  by  a  Saval  24  volt  D.  C.  solenoid. 

A  photograph  of  the  driver  section  appears  in  Fig.  3.  Two  different  types 
of  material  have  been  used  as  diaphragms  for  the  weakest  shock  waves. 
First,  Dupont  220  MD-31  cellophane  was  used.  This  material  tended  to 
shatter  into  small  pieces  which  necessitated  frequent  swabbing  of  the 
shock  tube.  The  second  type  diaphragm  material  that  was  used  was  Dupont 
mylar,  .  0015  inches  thick.  The  mylar  was  a  tougher  material  and  was 
used  for  the  stronger  shock  runs.  Unlike  the  cellophane,  the  mylar  did 
not  shatter  when  punctured  but  just  folded  back.  However,  more  energy 
was  lost  in  this  folding  process  and  it  took  a  higher  driver  pressure  using 
mylar  to  achieve  the  same  shock  Mach  number  than  it  did  when  cellophane 
was  used. 
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Figure  3.  Shock  Tube  Driver  Section 


3.  2-2  Optical  Equipment 

The  test  section  windows,  through  which  the  x  versus  t  history  of  the 
particles  is  recorded,  are  11/2  inches  high  and  4  inches  long  of  optically- 
flat  glass. 

A  modified  AVCO  rotating  drum  camera  is  used  in  conjunction  with  an 
oscillating  light  source  to  take  shadowgraph  pictures  of  the  particle's 

trajectory.  The  physical  layout  of  the  optical  system  is  reproduced  in 
Fig.  4. 

The  drum  camera  is  mounted  on  a  lathe  bed  with  two  compound  rests, 
with  rotating  drum  and  main  camera  body  on  one  compound  rest  and  the 
camera  lens  mounted  on  the  other.  A  flexible  bellows  connects  the  lens 
and  the  camera  body.  This  arrangement  permits  accurate  and  independent 
movement  of  the  lens  and  body  both  perpendicular  and  parallel  to  the  test 
section.  Thus  the  magnification  can  be  changed  as  desired  and  the  length 
of  the  test  section  can  be  traversed  with  the  lathe  setup.  The  lens  is  a 
Goerz  Red  Dot  Artar  process  lens,  which  is  a  highly  corrected  lens  de¬ 
signed  for  applications  where  the  image  to  film  distance  and  the  object  to 
film  distance  are  of  the  same  order. 

The  light  source  is  a  PEK  XE-9  Xenon  flash  tube.  The  energy  is 
switched  to  the  Xenon  tube  via  a  modified  version  of  a  Edgerton, 
Germeshausen,  and  Grier  LS-10  multiple  microflash  system.  This  system 
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Figure  4.  Optical  System 
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includes  a  low  voltage  power  supply,  a  pulse  shaper  unit,  a  time  delay 
unit,  and  five  high  voltage  discharge  units.  The  switching  circuits  con¬ 
trol  the  discharge  to  the  Xenon  flash  tube  of  one  .  001  microfarad  capacitor 
bank  and  four  .  005  microfarad  capacitor  banks,  all  charged  to  12  kilovolts. 
The  switching  circuit  is  activated  by  a  signal  generator  pulse-shaper  unit. 
By  varying  the  frequency  of  the  signal  generator,  it  is  possible  to  vary  the 
time  interval  between  flashes  from  10  milliseconds  to  10  microseconds. 

A  schematic  diagram  of  the  modified  EG  and  G  unit  is  in  Fig.  5. 

To  achieve  the  shortest  possible  rise  time,  it  is  necessary  to  reduce 
the  circuit  inductance  to  a  minimum  since  the  resistance  of  the  ionized 
Xenon  gas  is  only  2  fl .  This  was  accomplished  by  modifying  the  high 
voltage  discharge  units  of  the  multiple  microflash  system.  The  high 
voltage  section  of  the  first  switching  unit  was  installed  as  close  to  the 
Xenon  tube  as  possible  as  is  in  Fig,  6.  A  spark  gap  is  used  in  the  first 
unit  as  the  switch  for  the  energy.  This  was  done  since  the  spark  gap  will 
switch  the  current  faster  than  the  mercury  diode  of  the  other  four  units. 
With  the  above  arrangement  it  was  possible  to  obtain  a  480  nanosecond 
width  for  the  light  pulse  to  decay  to  10  per  cent  of  its  peak  value. 

The  Xenon  flash  lamp  is  mounted  on  a  drill  press  milling  vise,  per¬ 
mitting  motion  bcth  perpendicularly  and  parallel  to  the  test  section.  Both 
collimated  and  quasi -collimated  light  have  been  used  in  the  system. 
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Figure  5.  Schematic  Diagram  of  Microflash  System 
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First  high  voltage  Xenon  flash  tube 

switching  unit 

Figure  6.  First  Switching  Unit  and  Zenon  Flash  Tube 
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Two  types  of  film  have  been  used  in  the  experiment:  Kodak  Royal  Pan 
and  Kodak  Plus  X  film. 

3.2-3  Particle  Injector  System 

Within  the  particle  injector,  shown  in  Fig.  7,  the  particles  are  placed 
on  a  circular  platform  which  is  rotated  on  its  axis  by  means  of  a  solenoid 
thus  injecting  the  particles  into  the  test  section.  The  particle  injector 
is  designed  so  that  its  pressure  is  the  same  as  that  in  the  test  section  and 
the  hole  in  the  test  section  is  .  059  inches  in  diameter. 


Particle  injector  Solenoid 


Figure  7.  Particle  Injector 


3t  2-4  Shock  Speed  Measurement 

The  shock  speed  is  measured  with  two  flush-mounted  Kistler  model 
701A  quartz  pressure  transducers,  1. 1667  feet  apart,  in  conjunction  with 
a  microsecond  timer.  The  transducers  have  a  rise  time  of  ten  micro¬ 
seconds,  Two  Kistler  model  566  multi-range  electrostatic  charge 
amplifiers  are  used  to  amplify  the  signals  from  the  transducers.  These 
signals  are  in  turn  fed  to  the  start  and  stop  channels  of  a  Transistor 
Specialties  Incorporated  microsecond  Timer.  The  above  transducers 
are  very  sensitive  and  thus  allow  the  recording  of  very  weak  shock  waves. 
Desiring  the  best  possible  measurement  of  the  wave  speed,  the  transducers 
were  mounted  equally  upstream  and  downstream  of  the  test  section. 

3.  2-5  Pressurization  and  Vacuum  System 

A  schematic  diagram  of  the  pressurization  and  vacuum  system  is 
shown  in  Fig.  8.  The  system  was  designed  so  that  the  shock  tube  and 
driver  section  could  both  be  evacuated  to  any  desired  pressure.  The 
valve  arrangement  is  such  that  the  initial  pressure  can  be  controlled  in 
the  driver  section  and  shock  tube  independently.  In  a  similar  fashion 
both  the  driver  section  and  the  shock  tube  can  be  pressurized  in  an  inde¬ 
pendent  controlled  fashion.  With  this  flexibility  it  is  possible  to  maintain 
certain  desired  shock  strengths  while  varying  the  Re  number  and  M  /  VRe. 

X\ 
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Figure  8.  Schematic  Diagram  of  the  Pressurization  and  Vacuum  System 
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3.  2-6  Sequential  Timing  During  Experiment 

To  have  the  particles  in  correct  position  and  to  begin  the  photographing 
sequence  at  that  time,  an  accurate  timing  sequence  is  necessary.  Figure  9 
is  a  schematic  diagram  of  this  sequencing  system.  Once  the  particles  are 
dropped,  there  must  be  an  accurate  delay  before  the  solenoidal  drive  rod 
punctures  the  driver  diaphragm.  This  is  accomplished  by  means  of  a  thyra- 
tron  time  delay  unit.  A  circuit  diagram  of  this  circuit  appears  in  Fig.  10. 

A  second  sequencing  unit  is  needed  to  start  the  photographic  process. 

This  is  accomplished  by  having  the  signal  from  the  first  pressure  trans¬ 
ducer  sent  to  a  second  time  delay  unit  which,  after  a  sufficient  delay  al¬ 
lowing  time  for  the  shock  to  travel  to  and  interact  with  the  particles,  sends 
another  signal  to  the  first  high  voltage  switching  circuit.  The  high  voltage 
switching  circuit  then  dumps  its  energy  into  the  XE-9  flash  tube  and  the 
picture -taking  sequence  is  begun.  Figure  11  is  a  view  of  the  overall 
facility. 

3.2-7  Experimental  Procedure 

The  normal  experimental  procedure  is  as  follows: 

1.  Load  camera 

2.  Record  atmospheric  pressure 

3.  Record  atmospheric  temperature 

4.  Install  diaphragm 
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Figure  9.  Schematic  Diagram  of  Time  Delay  Sequencing  System 
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Figure  10.  Circuit  Diagram  of  Time  Delay  Unit 


Figure  11.  General  Shock  Tube  Facility 


5.  Load  injector 

6.  Evacuate,  pressurize,  or  both 

7.  Adjust  sensitivity  levels 

8.  Read  test  section  pressure 

9.  Run  rotating  drum  camera  up  to  speed 

10.  Engage  start  switch 

11.  Record  microsecond  timer  reading 

12.  Mark  and  remove  film  from  camera 

13.  Develop  film. 

3.  3  CALIBRATION  OF  EQUIPMENT 

All  of  the  parameters  except  p^  which  are  needed  to  calculate  are 
obtained  from  data  that  is  recorded  during  an  experiment.  Thus  it  is  es¬ 
sential  that  the  experimental  equipment  has  been  properly  calibrated,  and 
this  calibration  procedure  is  described  below. 

3.  3-1  Shock  Velocity  Measurement 

The  convective  flow  velocity,  Ug,  and  the  gas  density  in  the  convective 

flow  region,  p are  obtained  from  the  measurements  of  shock  velocity, 

C  ,  initial  temperature,  and  initial  pressure  and  application  of  the  normal 
s 

shock  relations.  For  weak  shock  waves  is  very  sensitive  to  changes  in 

C  ,  or  the  time  for  the  shock  wave  to  travel  the  distance  between  the  two 
s 

pressure  pickups.  This  time  measurement  must  therefore  be  very  accurate. 
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The  linearity  and  slope  of  the  signals  from  the  pressure  transducers 
were  checked  by  putting  the  output  from  the  amplifiers  of  the  transducers 
into  a  Tektronix  oscilloscope  and  photographing  the  trace  with  a  Polaroid 
attachment.  Photographs  of  the  traces  were  taken  of  the  start  transducer 
and  its  amplifier  and  of  the  stop  transducer  and  its  amplifier  at  equal 
shock  strengths.  The  amplifiers  were  then  interchanged  and  the  procedure 
was  repeated  at  the  same  shock  strengths.  Upon  comparing  the  results 
it  was  found  that  the  four  traces  are  virtually  indistinquishable.  Two 
pictures  of  such  traces  appear  in  Fig.  12a  and  12b. 

The  sensitivity  adjustments  on  the  start  and  stop  channels  of  the 
microsecond  timer  were  calibrated  so  that  both  channels  could  be  set  at 
the  same  level.  The  gain  on  the  amplifier  and  the  sensitivity  on  the 
microsecond  timer  were  always  adjusted  so  that  both  channels  of  the 
timer  would  be  activated  on  about  the  first  10  per  cent  of  the  signal  from 
the  pressure  transducers. 

Since  the  slopes  and  linearity  of  the  two  transducer  amplifier  com¬ 
binations  were  the  same  and  the  sensitivities  on  the  timer  were  adjusted 
to  the  same  level,  the  shock  speed  that  was  measured  by  the  transducer 
microsecond  timer  combination  would  be  the  actual  shock  speed  if  shock 
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speed  decay  was  unimportant.  By  setting  the  same  driver  to  test  section 
pressure  ratio  and  making  a  series  of  shock  runs,  the  time  interval  read 
on  the  microsecond  timer  would  consistently  be  within  plus  or  minus  one 
microsecond  out  of  900  microseconds. 

3.  3-2  Schlieren  Photographs  of  Shock  Front 

A  Schlieren  system  was  set  up  to  verify  the  existence  of  a  plane  shock 
front  and  to  make  sure  that  the  injector  hole  did  not  cause  the  front  to 
bend  locally.  A  spark  gap  triggered  by  the  first  unit  of  the  multi-microflash 
was  used  as  the  light  source  for  the  Schlieren  system.  By  incorporating 
the  spark  gap  with  the  microflash  system  and  varying  the  time  delay  in 
the  microflash  system,  it  was  possible  to  take  photographs  of  the  shock 
front  at  various  locations  along  the  viewing  area  of  the  test  section.  In  all 
cases  the  shock  front  was  plane  and  perpendicular  to  the  direction  of  flow. 
The  fact  that  the  shock  front  is  plane  and  is  perpendicular  to  the  flow  di¬ 
rection  is  justification  for  using  the  normal  shock-tube  relationships. 

3.  3-3  Calibration  of  Optical  Equipment 

The  calibration  of  the  optical  equipment  is  extremely  important  for 
both  the  particle  diameter,  and  the  distance  versus  time  data  are  obtained 
photographically.  The  camera  was  focused  and  the  magnification  was 
determined  in  the  following  way.  A  cylindrical  rod,  .  043  inches  in 
diameter,  was  inserted  through  the  particle  injector  hole  parallel  to  the 
test  section  walls.  The  rotating  drum  camera  was  then  focused  on  the  rod, 
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and  a  picture  of  the  rod  was  taken  with  the  drum  stationary.  The  film  was 
developed  and  the  negative  was  then  viewed  through  a  Bausch  and  Lomb 
microscope  with  graduated  eyepiece  to  determine  the  size  of  the  image. 

The  magnification  of  the  microscope  had  previously  been  determined  using 
the  grid  from  an  Edmund  comparator.  Thus  the  magnification  of  the  camera 
was  known.  Some  doubt  remained  as  to  whether  this  magnification  was 
correct  in  that  the  calibration  probe  was  six  times  larger  in  diameter 
than  the  spherical  particles  and  the  test  object  was  not  spherical  like  the 
particles.  Accordingly,  synthetic  sapphire  spheres,  .015  inches  in  dia¬ 
meter  (±  .  0001  inches),  were  dropped  through  the  particle  injector  and 
photographed.  The  photographic  images  of  the  sapphire  balls  yielded  sphere 
diameters  which  were  within  .  1  per  cent  of  the  actual  sapphire  sphere 
diameter.  Thus,  the  method  for  determining  the  magnification  factor 
proved  acceptable. 

Another  optical  test  was  made  to  determine  whether  the  lines  on  a 
grid  pattern  remained  straight  and  undistorted  so  that  the  particle  displace¬ 
ment  measurements  would  be  true.  A  grid  network  on  a  clear  plastic 
was  inserted  into  the  test  section  and  photographed.  By  viewing  the  re¬ 
sulting  image  on  the  film  no  distortion  could  be  detected.  Thus  the  dis¬ 
placement  measurements  on  the  film  would  be  a  true  representation  of 
the  actual  movement  of  the  particle. 
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To  determine  the  influence  of  development  time  on  image  size,  a  final 
optical  test  was  made.  Three  wires,  0.  0036  inches,  0.  0070  inches,  and 
0. 0104  inches  in  diameter,  were  photographed  nine  times  under  the  same 
conditions.  Using  Kodak  D-ll  developer  at  68°F,  each  piece  of  exposed 
film  was  developed  at  times  ranging  from  five  to  nine  minutes  at  half- 
minute  invervals  —  seven  minutes  being  normal.  Upon  microscopic 
examination,  no  apparent  diameter  change  could  be  detected  in  any  of  the 
three  wires.  Thus  in  actual  experimental  runs,  it  was  determined  that 
the  development  time  could  be  varied  plus  and  minus  one  minute  without 
affecting  the  validity  of  the  experimental  results. 

3.  3-4  Shock  Tube  Attenuation 

Ideally  a  shock  wave  propagates  at  a  constant  velocity  in  a  shock  tube. 

In  actuality,  however,  the  shock  wave  attenuates.  The  amount  of  this 
attenuation  depends  on  several  things:  the  strength  of  the  original  wave 
at  the  diaphragm,  the  distance  of  the  measuring  device  from  the  dia¬ 
phragm,  and  the  size  of  the  shock  tube.  An  experimental  study  of  at- 

(15) 

tenuation  in  shock  tubes  was  conducted  by  R.  J.  Emrich  and  C.W.  Curtis' 
They  conducted  their  experiments  in  shock  tubes  of  various  sizes.  For  a 
shock  tube  with  the  same  hydraulic  radius  as  in  the  present  study,  they 
found  the  following  relationship  to  be  valid  within  a  factor  of  two. 


dM 

dx 


-  5.  08  x  10 
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The  shock  Mach  number  was  always  less  than  1.  25.  For  M  =  1.25  dM/dx 


is 


For  x  =  14  inches 


AM 

Ax 


-2.  3  x  10 


AM  -  .  003 


Thus  the  attenuation  over  the  interval  in  which  the  velocity  is  measured 
is  negligible  and  the  convective  flow  conditions  can  be  calculated  directly 
without  any  correction  for  attenuation. 

3.  3-5  Shock  Tube  Test  Time 

The  shock  tube  test  time  depends  on  many  factors:  driver  section 
length,  the  distance  between  the  driver  section  and  the  test  section,  the 
distance  between  the  test  section,  and  end  of  the  shock  tube,  the  particular 
gases  or  gas  being  used,  the  strength  of  the  shock,  and  the  local  speed  of 
sound.  The  test  time  can  be  determined  analytically,  however,  the  actual 
test  time  is  usually  less  and  thus  a  series  of  experiments  were  performed 
at  various  shock  strengths  to  determine  the  actual  test  time.  This  was 
accomplished  by  putting  the  output  of  the  pressure  pickups  into  a 
Tektronix  oscilloscope  and  photographing  the  traces.  The  first  pressure 
pickup  started  the  sweep  of  the  oscilloscope.  A  typical  trace  is  shown  in 
Fig.  13.  The  test  time  extends  from  the  instant  of  pressure  increase  to 
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Figure  13.  Shock  Tube  Test  Time 


the  onset  of  the  rarefaction  wave  and  from  the  above  figure  is  4.  4  milli¬ 
seconds.  The  second  pressure  pulse  is  from  the  reflected  shock  wave. 
The  shortest  run  time  that  was  recorded  over  the  desired  range  of  shock 
strengths  was  4.  2  milliseconds.  Thus  all  of  the  experimental  runs  were 
made  in  4.2  milliseconds  or  less. 

i 

i 

| 
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IV.  DATA  REDUCTION  AND  ERROR  ANALYSIS 


4. 1  RANGE  OF  EXPERIMENTAL  MEASUREMENTS 

The  experimental  study  consisted  of  determining  the  drag  coef¬ 
ficient  of  non-burning  particles  under  unsteady  flow  conditions.  For  all 
the  experimental  runs,  the  relative  Mach  number  of  the  flow  was  always 
less  than  .  3.  The  M^/  VRe  was  such  that  the  flow  was  always  in  the  con¬ 
tinuum  regime.  The  Reynolds  number  ranged  from  150  to  1700. 

Three  different  types  of  particles  were  used  in  the  study:  glass 
beads,  Winchester  Western  HP  295  ball  powder,  and  synthetic  sapphire 
balls.  The  superbrite  glass  beads  were  obtained  from  Reflective  Products 
Division  of  the  Minnesota  Mining  and  Manufacturing  Company.  The  glass 
beads  for  the  most  part  were  spherical;  however,  some  discretion  was 
needed  in  examining  the  photographs  of  the  glass  beads  to  sort  out  the 
few  beads  which  were  not  spherical.  The  Winchester  Western  HP  295 
ball  powder  likewise  was  mostly  spherical,  but  again,  some  selection  of 
data  was  necessary  to  eliminate  the  non-spherical  powder.  The  preci¬ 
sion-lapped  sapphire  balls,  obtained  from  Industrial  Tectonics,  Inc. , 
were  very  spherical  (within  0.  000010).  The  balls  are  also  uniform  in 
their  diameter  (±  0.0001  inches),  and  their  surface  finish  is  precise, 

1.  5  microinches.  A  table  listing  some  of  the  particle  characteristics 
appears  below: 
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TABLE  I 


Type  of  Particle 

Density 

(gm/cc) 

Size  Range  or  Exact  Size 

Superbrite  Glass 

2.49 

150  ix  -  250  i a 

Beads 

Winchester -We  stern 

1.67 

280  ix  -  350  ix 

HP  295  Ball  Powder 

Sapphire  Balls 

3.978 

396. 8  ix 

The  density  of  the  Superbrite  glass  beads  were  obtained  experimentally 
using  a  Beckman  Air  Comparison  Pycnometer. 

4.  2  REDUCTION  OF  EXPERIMENTAL  DATA 

From  each  experimental  run  the  following  variables  are  needed  for 
the  determination  of  C^  and  Re 

1.  particle  diameter 

2.  convective  flow  velocity 

3.  convective  gas  density 

4.  gas  viscosity  in  the  convective  flow  regime 

5.  particle  velocity 

6.  particle  acceleration. 

As  mentioned  previously,  the  particle  diameter  is  obtained  by  viewing 
its  image  on  the  film  under  a  microscope.  Knowing  the  camera  and  micro¬ 
scope  magnification  factors,  the  particle  diameter  is  obtained  directly. 
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The  convective  flow  velocity  and  gas  density  are  obtained  from  the  normal 
shock  relations  knowing  the  shock  speed,  the  temperature  of  the  undis¬ 
turbed  gas,  and  the  pressure  of  the  undisturbed  gas  in  the  test  section. 

The  temperature  in  the  convective  region  is  likewise  given  by  the  normal 
shock  relations.  The  viscosity  values  for  air  at  various  temperatures 
were  obtained  from  air  viscosity  tables  of  the  National  Bureau  of  Standards 
For  small  temperature  ranges,  i.  e. ,  AT  ~  25°R,  the  viscosity  versus 
temperature  curve  is  linear.  Thus,  a  series  of  curves  of  viscosity  versus  tem¬ 
perature  were  plotted  for  different  temperature  ranges  from  the  above 
table.  Linear  viscosity  relationships  versus  temperature  were  obtained 
from  these  curves  for  each  temperature  range  so  that  the  viscosity  could 
be  solved  for  analytically.  The  remaining  two  variables  which  must  be 
found  for  each  particle  are  the  velocity  and  acceleration. 

The  particle  distance  versus  time  history  is  recorded  by  taking  five 
photographs  on  the  rotating  drum  camera  at  equal  time  intervals.  From 
these  photographs,  the  particle  velocity  and  acceleration  must  be  obtained. 

4.  3  METHOD  OF  DATA  REDUCTION 

The  simplest  technique  for  the  reduction  of  the  x  versus  t  data  would 
be  to  set  up  a  difference  table.  If  this  method  was  used,  one  would  obtain 
the  velocity  of  the  particle  from  the  first  difference  and  the  particle's 
acceleration  from  the  second  difference.  The  main  disadvantage  of  the 
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above  method  is  that  each  second  difference  depends  on  three  position 
measurements,  and  is  very  sensitive  to  any  errors  in  these  measure¬ 
ments.  Another  disadvantage  is  that  an  average  acceleration  and  velocity 
is  obtained  rather  than  an  instantaneous  value.  It  would  be  more  accurate 
to  fit  the  position  data  with  a  polynominal.  Since  there  are  five  position: 
points,  the  highest  order  polynominal  which  can  be  fitted  to  the  points 
is  a  fourth  order  polynominal  which  would  pass  through  each  position 
point.  Thus,  any  measurement  error  in  the  position  data  would  be  ampli¬ 
fied  when  the  acceleration  was  obtained  by  differentiating  the  analytic  x 
versus  t  expression  twice,  since  a  fourth  order  polynomial  would  pass 
through  every  point. 

A  second  order  polynomial  could  be  used  to  fit  the  position  data. 
However,  this  implies  that  the  force  on  the  particle  is  a  constant.  The 
velocity  relative  to  the  particle  is  actually  changing  with  time  and  thus  the 
force  on  the  particle  is  changing  with  time.  The  remaining  choice  is  a  third 
order  polynominal.  By  using  the  third  order  fit,  the  final  curve  would 
not  be  forced  to  pass  exactly  through  the  position  data  and  thus  some  mea¬ 
surement  error  could  exist  without  affecting  the  final  result  drastically. 

A  least  square  method  was  used  to  fit  a  third  order  curve  through  the 
five  position  points.  Once  the  analytic  expression  was  obtained,  the  velo¬ 
city  and  acceleration  of  the  particle  with  respect  to  inertial  space  was 
obtained  by  differentiating  the  x  versus  t  expression.  The  time  which  was 
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substituted  in  the  analytic  expression  for  x  was  that  which  corresponded 
to  the  third  picture.  It  was  felt  that  the  second  derivative  was  more  cor¬ 
rect  near  the  center  point  of  the  x  versus  t  expression  than  at  the  times 
corresponding  to  the  first  or  fifth  pictures. 

The  data  reduction  process  for  the  most  part  is  handled  by  an  IBM  7090 
computer  program.  The  data  which  is  put  into  the  computer  program  for 
each  run  is  as  follows: 

1.  Run  number 

2.  Particle  number 

3.  Diameter  of  particle  on  film 

4.  Position  data;  Xq,  x^,  x^,  x^,  x^ 

5.  Atmospheric  pressure 

6.  Initial  temperature 

7.  Initial  test  section  pressure 

8.  Oscillator  frequency 

9.  Microsecond  timer  reading 

10.  Magnification 

11.  Density  of  particles 

With  the  above  data  the  computer  calculates  p^,  Tg,  Ug,  and  solves 
for  the  coefficients  of  the  third  order  polynomial;  then  differentiates  this 


expression  to  obtain  the  particle  velocity  and  acceleration.  Finally,  it 
solves  for  and  Re  and  prints  this  and  other  pertinent  data  on  the 
output  page. 

A  second  data  reduction  technique  which  calculates  a  mean  was 
used  to  check  the  third  order  method  results  as  described  below.  The 
particle  acceleration  and  velocity  can  be  written  as 


a  = 


•  » 

=  x 


V 

P 


The  equation  for  can  be  written  as 

B(U2  -  x)2  =x 


(8) 


where 


B  = 


^2CD 

V 


Let  x  =  p  x  =  dp/  dt.  Equation  (8)  becomes 


dj? 

<U2  -  p)2 


Bdt 


(9) 


If  we  assume  that  over  the  time  interval  of  interest  B  is  a  constant  we  may 
integrate  Eq.  (9).  The  values  of  p  and  x  at  limits  of  the  integral  are 


t  =  0 
t  =  t 


P  =  P, 


P  =  P 
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Equation  (9)  becomes 


P  d(U2  -  p)  * 

I  - - - 2  =  B  I  dt 

pJ  <U2  -  oJ 


Performing  the  integration  and  substituting  the  limits 


P  =  U 


2  Bt  +  1 


dx 

dt 


U2‘po 


(10) 


The  above  may  be  rewritten  as 


x 


/  dx=f  u2dt-i  J t 


X 


o 


o 


B/U 


-pU  +  l 


J  ^B|U2-po)t  +  1l 


Integrating  and  substituting  in  for  the  limits  one  has 


x  “  x0  +  V  -  B  ln 


B(U2  -  VD  >  ‘  +  1 
po 


(11) 


A  least  squares  technique  is  used  on  the  position  data  to  solve  for  the 

coefficients  of  Eq.  (11),  i.  e. ,  x  ,  B,  and  V  .  Once  these  have  been 

°  P° 

evaluated,  can  be  obtained  from  B.  The  particle  velocity,  V  ,  was 
calculated  from  Eq.  (10)  using  the  time  corresponding  to  the  third  picture. 
Knowing  V  ,  Ug,  d,  p 2,  and  lhe  Reynolds  number,  Re2,  can  be  cal¬ 
culated.  These  calculations  have  also  been  programmed  for  the  computer. 
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The  Cp  method  assumes  that  the  Cj^  is  constant  over  the  time  interval 
of  the  test,  an  approximately  true  assumption  in  the  actual  case  since  for 
a  typical  run  ACD/CD  q  1. 1  x  10-2  over  the  test -time  interval. 

4.  4  TYPICAL  EXAMPLE  OF  DATA  REDUCTION 

Figure  14  shows  five  sequential  pictures  of  a  typical  run  (Run  50E) 
where  the  particles  are  being  accelerated  by  the  convective  flow  behind 
a  shock.  The  x  versus  t  measurements  were  made  using  a  ruler  graduated 
in  0. 01  of  an  inch  in  conjunction  with  an  Edmund  12-power  comparator. 
With  this  system  it  was  possible  to  estimate  readings  on  the  scale  to 
0.  003  of  an  inch. 


Various  information  obtained  from  Run  50E  for  the  data  reduction 
is  listed  below : 

U  =  3.10 

t  =  6. 812  x  10  4sec 
P 

d  =  5. 166  x  10“4ft 


T  ,  =  76.  40  F 

atm 

Patm  =  29' 10  inches  H& 


t 

m 


=  18.  28  inches  Hg 
=  945.  0  x  10  **  sec 


The  above  data  when  fed  into  the  two  computer  programs  gave  the  following 
results 
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Figure  14.  Typical  Example  of  x  Versus  t  Data  of  a  Particle 
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Mx  =  1.088 
U2  =  159. 14  ft/ sec 
P2  =  •  0519  lbs/ ft3 
=  1135.  0  ft/  sec 
a2  =  1166.  9  ft/  sec 
P2  =  1568.8  lbs/ ft2 
P2  =  1.287  x  10  3  lb/ ft,  sec 

From  this  point  on  the  method  of  data  reduction  by  the  computer  pro¬ 
grams  differed.  As  mentioned  previously,  the  first  program  used  a  third 
order  fit  to  the  x  versus  t  data.  A  graph  of  the  third  order  fit  for  Run  50E 
appears  in  Fig.  15.  The  particle^  velocity  and  acceleration  are  obtained 
by  differentiating  the  analytic  expression  for  x  versus  t  once  and  twice 
respectively.  This  expression  for  the  particular  particle  in  question  is 

x  —  4.  605  t  +  4.  223  x  103  t2  -  7.  045  x  104  t3 

The  time  which  was  substituted  in  for  the  velocity  and  acceleration  is 
that  which  corresponds  to  the  third  picture  and  is 

-3 

t2  =  1.  36  x  10  seconds 

The  velocity  and  acceleration  corresponding  to  t2  is 

V  =15.  72  ft/ sec 
P 

a  =  7869.  35  ft/ sec2 
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rANCE  IN  FE 


When  all  of  the  above  values  are  put  into  the  equations  for  CD  and  Re  one 
has 

CD  =  .  7926 
Re  =  298.81 

The  second  computer  program  solves  by  the  method  of  least  squares  for 

the  unknown  constants  x  ,  B,  and  V  in  Eq.  (11).  Once  these  quantities 
_  0  po 

are  known  C_  and  Re0  are  determined.  For  this  sample  run  one  gets 
D  Ci 

CD  =  .  7892 
Re2  =  299.  39 


In  the  present  use,  and  in  the  evaluation  of  other  data,  the  two 
methods  of  data  reduction  were  found  to  be  in  substantial  agreement. 

4.  5  GENERAL  EQUATIONS  FOR  DETERMINING  PROBABLE  ERROR 
In  order  to  assess  the  validity  of  the  experimental  data  an  error  ana¬ 


lysis  must  be  performed.  As  can  be  seen  from  Eq.  (5),  is  a  function 

of  the  following  independent  variables;  Pp,  a,  d,  p2,  Ug,  and  V^.  The 

(17) 

probable  error  in  can  be  written  as 


c  \2  ,  c  \2 

P<CD>=  \TT  (6“>  + 

\  P  /  '  ' 


ecf  ,  isc  \2  „ 

srj  <5d>  +yw 


jSC  \2  ,  C  f  2 

+  auf  (6U2)2+K  <sV 

\  2/  \  P 
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Performing  the  necessary  differentiation  and  substituting  into  Eq.  (12)  the 


probable  error  in  becomes 


P(CD) 

CD 


vv  i 


6_a 
a  i 


6p. 


+  4 


6U, 


U 


R 


'6  V 


+  4 


U 


R 


The  probable  deviation  or  errors  in  p  ,  a ,  d,  and  V  can  all  be  estimated 

P  P 

or  obtained  directly  from  the  data.  However,  the  errors  in  U2  and  p2 

must  be  evaluated  from  the  normal  shock  relationships  in  conjunction 

with  C  ,  a. ,  and  P. ,  i.  e. , 
s  i’  r  ’ 


u 


2 


rs _ 

(r  +  i)  cg 


(14) 


p2 


r  (r  +  i) 
f(r  -  l)  Cs2 


(15) 


In  all  the  calculations  y  =  1.4,  6  U2  and  6  p2  are  obtained  by  taking  dif¬ 
ferentials  of  Eqs.  (14)  and  (15).  The  actual  values  of  6  U2  and  6  p^  for  a 
particular  run  are  then  obtained  by  substituting  in  the  numerical  values 
for  Cg,  a^,  and  P^.  These  values  can  be  put  into  Eq.  (13)  with  the  other 
deviations  and  the  probable  error  in  drag  coefficient,  P(C^),  for  a  parti¬ 
cular  case  can  be  obtained. 
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After  differentiating  Eq.  (6)  and  substituting  the  results  into  Eq.  (16) 
P(Re),  becomes 


As  in  the  case  of  C^,the  deviations  6  V  ,  6  d,  and  6  p  can  be  estimated  or 

D  P  ^ 

obtained  directly  from  the  data.  The  deviations  6  p ^  and  6  Ug  have  been 
determined.  Thus  for  a  particular  run  the  probable  error  in  Reynolds 
number  can  be  obtained  using  Eq.  (17). 

4.  6  TYPICAL  EXAMPLE  OF  PROBABLE  ERROR 

The  same  run  which  was  used  in  the  data  reduction  section  as  an 
example  run  will  be  used  in  determining  the  probable  error  in  CD  and  Re, 
namely  Run  50E.  The  necessary  data  from  Run  50E  is  given  below: 
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Pp  =  2 .  475  gm/  cc 

t 

a  =  7869.  3  ft/ sec' 


d  =  5. 166  x  10  ft 


p2  =  .  0519  lb/ ft 

U2  -  159. 14  ft/ sec 

UD  =  143.43  ft/ sec 
K 


V  =15.72  ft/ sec 
P 

C  =  1234.  5  ft/  sec 
s 

=  1135.  0  ft/  sec 

=  1292.8  lb/ft2 
-fi 

Time  =  945.  0  x  10  sec 

P2  =  1. 2868:  x  10  ^  lb/  ft  sec 


The  probable  error  in  will  be  computed  first.  The  first  term  in  Eq.  (13) 

is  6  p  /  p  .  The  volume  of  the  glass  beads  was  found  by  using  a  Beckman 
P  P 

model  930  air  comparison  pycnometer.  A  10  c.  c.  or  larger  sample  of 
beads  is  needed  for  an  accurate  volume  determination.  Once  the  volume 


was  determined  the  sample  was  weighed  on  a  beam  balance  in  order  to  ob¬ 
tain  the  average  density.  The  6  p^  and  average  density,  p^  using  the  above 
technique  was 


from  which 


«P. 


=  .  050  gm/  cc 
=  2.  475  gm/cc 

=  4.  08  x  10"4 


Consider  the  second  term  in  Eq.  (13),  6  a/ a .  The  probable  error  in 
determining  acceleration  was  attributed  only  to  incorrect  particle  dis¬ 
placement  measurements.  Typical  displacement  errors  (±  .  003  inches) 
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were  incorporated  into  the  original  measurements  of  x  versus  t.  This 
new  data  was  then  run  through  the  computer  program  for  data  reduction 
to  find  the  new  value  for  acceleration.  Doing  this  the  deviation  in  ac¬ 
celeration  for  Run  50E  is 

6  a  =  274.  5  ft/  sec^ 
j^j  =  12. 18  x  l(f 4 

The  next  term  is  6  d/ d.  The  deviation  in  particle  diameter  was  con¬ 
sidered  due  only  to  the  error  in  estimating  the  exact  location  of  the 
particle's  edge.  The  deviation  for  the  diameter  was 

6  d  =  .  151  x  10  4  ft 


=  8.  52  x  10 


The  deviation  in  the  particle's  velocity,  6  V  ,  was  calculated  as  being 

P 

due  to  incorrect  displacement  measurements  as  in  the  acceleration  case. 

The  method  for  determining  6  was  entirely  analogous  to  that  of  5  a  as 

mentioned  above,  and  for  Run  50E  6  V  is 
’  P 

5  V  =  .24  ft/ sec 
P 


11  x  10 
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6  pg  depends  on  the  deviations  of  Cg,  a^,  and  P^.  The  initial  test  section 
pressure,  P^,  is  read  off  a  mercury  manometer  and  thus  6  P^  is  rather 
small. 


6PX  =  5.17  lb/ ft 


6P 


-=  .  4x  10‘2 


The  deviation  in  a^  is  that  caused  by  a  change  in  the  temperature  and  is 
for  Run  50E 

6  a^  =  2  ft/  sec 

The  deviation  in  shock  speed,  C  ,  is  considered  due  only  to  the  error  in 

s 

the  microsecond  timer.  It  was  shown  in  a  previous  section  that  attenua¬ 
tion  affects  could  be  neglected.  It  was  also  mentioned  that  since  the  outputs 
of  the  two  transducers  had  the  same  slopes  and  both  channels  of  the  micro¬ 
second  timer  were  calibrated,  the  measured  time  would  be  correct.  The 
deviation  due  to  the  error  in  the  microsecond  timer  is 

6  C  =  1. 2  ft/  sec 
s 

Substituting  in  the  proper  values  in  6  p ^  we  have 

6  Pg  =  -  .  416  x  10  ^ 
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-  =  .  645  x  10-6 


Similarly  6  Ug  is  obtained  and  is 


6U2  =  -1.226  ft/ sec 


/6  U2  -4 

4j-rj— —  =2.92x10* 

\R 


All  of  the  terms  which  appear  in  Eq.  (14)  have  been  computed.  Putting 
these  terms  in  Eq.  (13)  the  probable  error  in  drag  coefficient  is  for 
Run  50E 


P(CD}  -2 

=  5.28  x  10 


In  calculating  the  probable  error  in  Reynolds  number  Eq.  (17)  must 
be  used.  This  equation  contains  the  same  deviations  as  those  just  cal¬ 
culated  for  C  plus  6  ii  .  Thus  for  the  Re  number  case,  it  is  necessary 
D  ^ 

only  to  find  6  ju 5  ii^  is  caused  only  by  an  error  in  the  temperature 
measurement.  This  gives  6  as 


6  /i2  =  .  005  x  10"  lb/ ft,  sec 


6  /z =  1.  51  x  10  ^ 
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Putting  (6  jUg/  j Ug)  and  the  other  values  into  Eq.  (17)  for  the  probable 
error  in  Re  for  Run  50E  we  have 


P(Re) 

Re 


=  3.71  x  10 
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V.  RESULTS  AND  DISCUSSION  OF  RESULTS 
5. 1  EXPERIMENTAL  RESULTS 

The  experiments  involving  the  glass  beads  were  all  carried  out  using 
a  magnification  of  3. 1.  Before  a  particular  particle  was  used,  the  film 
record  of  the  x  versus  t  history  was  checked  against  the  following  condi¬ 
tions. 

1.  Any  particle  which  was  closer  than  five  diameters  to  another 
particle  perpendicular  to  the  direction  of  motion  was  not  used. 

2.  Any  particle  which  was  within  twelve  diameters  of  another 
particle  wake  or  had  another  particle  within  twelve  diameters 
of  its  wake  was  not  used  in  the  data  reduction. 

3.  If  any  diaphragm  material  appeared  in  any  of  the  five  frames, 
the  run  was  discarded. 

4.  Only  spherical  particles  were  used  in  the  data  reduction. 

5.  Only  particles  which  were  sufficiently  sharp  were  used  for  the 
data  reduction. 

Because  of  the  above  restrictions,  many  particles  could  not  be  used  to 
obtain  CD  versus  Re  data.  Table  2,  in  the  Appendix,  is  a  summary  of 
experimental  data  for  glass  beads  in  incompressible  continuum  flow. 

CD  and  Re  results  obtained  using  both  methods  of  data  reduction  appear 
in  the  table.  If  the  versus  Re  points  are  plotted  on  log-log  paper, 
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they  appear  to  fall  in  nearly  a  straight  line.  Thus,  a  least  squares  technique 
was  used  to  fit  the  versus  Re  data  with  an  equation  of  the  following  form: 

2 

log10  Cd  ~  Ao  +  A1  log10  Re  +  *2  ^log10  Re^  (18) 

Other  analytical  forms  were  tried  for  versus  Re  data,  however,  none 
of  these  other  forms  fit  the  data  nearly  as  well  as  that  of  Eq.  (18).  The 
results  which  appear  in  Table  2  are  plotted  in  Figs.  16  and  17  for  the 
two  methods  of  data  reduction.  The  best  fit  curve  for  these  points  is  also 
plotted. 

A  second  series  of  experimental  runs  were  made  with  the  glass  beads 
at  higher  shock  strengths  to  begin  a  study  of  the  influence  of  compressi¬ 
bility.  The  results  of  these  runs  appear  in  Table  3.  The  relative  Mach 
number,  M  ,  for  the  flow  field  about  the  particles  is  .  15  <  M_,  <  .  30. 

The  has  increased  as  much  as  50  per  cent  with  respect  to  similar 
data  which  appeared  in  Table  2.  The  scatter  in  for  a  particular  Re 
has  increased  approximately  40  per  cent.  This  shift  in  cannot  be  at¬ 
tributed  to  compressibility  effects  since  the  relative  Mach  number,  M_,, 
is  still  quite  low.  Likewise  the  scatter  in  cannot  be  explained  by  dif¬ 
ferences  in  experimental  technique  since  the  results  which  appear  in  Table 
3  were  obtained  in  precisely  the  same  manner  as  the  data  which  appears 
in  Table  2.  Thus  two  more  series  of  experimental  runs  were  made  in  an 
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curve  for  incompressible  continuum  non 
spherical  glass  particles 


Coefficients  Versus  Reynolds  Number  for  Third  Order  Polynomial  Data 
Reduction  Technique 


curve  for  incompressible  continuum  non 
spherical  glass  particles 


[ti] 


Figure  17.  Drag  Coefficient  Versus  Reynolds  Number  for  Mean  Drag  Coefficient  Method 


attempt  to  explain  the  above  results.  They  were  made  using  the  same 
experimental  technique  that  was  used  for  the  data  of  Tables  2  and  3.  The 
only  change  was  the  types  of  particles  that  were  used,  Winchester-Western 
HP  295  ball  powder  and  sapphire  balls.  The  data  summaries  for  the 
HP  295  ball  powder  and  sapphire  balls  appears  in  Tables  4  and  5  res¬ 
pectively.  The  results  of  Tables  3,  4,  and  5  are  illustrated  in  Fig.  18. 

5.  2  SURFACE  ROUGHNESS 

Figure  18  shows  that  the  scatter  in  the  HP  295  ball  powder  data  is  of 
the  same  level  as  that  obtained  for  the  higher  relative  Mach  number  runs 
using  glass  beads.  The  relative  Mach  number  of  the  flow  about  the  ball 
powder,  however,  was  MR  <  .  124.  Both  for  the  HP  295  ball  powder 
and  the  higher  relative  Mach  number  glass  bead  data,  the  average 
for  a  particular  Reynolds  number  fell  above  the  best  fit  curve  of  the  glass 
beads.  In  contrast  to  this,  the  average  for  the  sapphire  runs  for  a 
particular  Reynolds  number  falls  below  the  best  fit  curve  of  glass  beads 

and  does  not  seem  to  be  a  function  of  the  MD. 

R 

As  an  attempt  to  explain  the  varied  results  which  appear  in  Fig.  18, 
photomicrographs  were  taken  of  the  particles.  Photomicrographs  of 
200  p.  glass  beads  appear  in  Figs.  19,  20,  and  21.  Under  the  relatively 
low  power  magnification  of  Fig.  19,  the  glass  beads  appear  to  have  fairly 
smooth  surfaces  except  for  a  few  bubbles  or  craters  which  appear  on 
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Best  fit  curve  for  incompressible 
continuum  non-burning  spherical 
glass  particles 


Figure  18.  Drag  Coefficient  Versus  Reynolds  Number — Glass  Particles,  Ball,  Powder, 

Sapphire  Balls 


Figure  19.  Photomicrograph  of  Glass  Particles,  Magnification  =  200 
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Figure  20.  Photomicrograph  of  Glass  Particle,  Magnification  =  1840 
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Figure  2h 


Photomicrograph  of  Glass  Particle,  Magnification  =  1840 
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most  of  the  beads.  Figure  20  shows  a  detailed  view  of  these  bubbles  or 
craters.  These  protrusions  and  craters  are  quite  pronounced  relative 
to  the  particle  diameter.  Figure  21  is  a  detailed  view  of  a  portion  of  the 
glass  bead  which  is  free  of  the  larger  surface  blemishes.  However,  this 
portion  of  the  surface  still  possesses  a  considerable  amount  of  rough¬ 
ness  relative  to  the  particle  diameter.  Photomicrographs  of  HP  295  ball 
powder  appear  in  Figs.  22,  23,  and  24.  Figure  22  is  a  low-power  photo¬ 
micrograph  and  shows  only  that  the  ball  powder  is  spherical  with  a  few 
surface  protrusions.  Figure  23  is  of  intermediate  magnification,  and  it 
appears  from  this  figure  that  the  entire  surface  of  the  ball  powder  is 
rough  relative  to  the  similar  photomicrograph  of  the  glass  beads  in 
Fig.  19.  Finally,  Fig.  24  is  a  detailed  view  of  the  ball  powder.  Again 
the  surface  appears  to  be  rougher  than  that  of  glass  beads.  Photomicro¬ 
graphs  of  sapphire  balls  appear  in  Fig.  25,  26,  and  27.  Figure  25  is  a 
low -power  photomicrograph  and  indicates  only  that  the  sapphire  ball  is 
spherical.  Figure  26  is  a  detailed  view  of  a  sapphire  ball.  From  this 
figure  it  is  evident  that  the  sapphire  balls  have  less  surface  roughness 
than  even  the  portions  of  glass  that  were  free  of  the  bubbles  or  craters. 
Figure  27  is  a  photomicrograph  of  a  sapphire  ball  which  was  washed  in 
trichloroethylene  before  the  photomicrograph  was  taken.  Almost  all  of 
the  roughness  which  appeared  in  Fig.  26  is  missing  in  Fig.  27.  The 

relative  Mach  number  range  for  the  sapphire  was,  .14  <  M  <  .313 

—  R  — 
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Figure  22.  Photomicrograph  of  HP  295  Ball  Powder,  Magnification  =110 
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Figure  23.  Photomicrograph  of  HP  295  Ball  Powder,  Magnification  =  190 
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Figure  25.  Photomicrograph  of  Sapphire  Balls,  Magnification  =  202 
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yet  none  of  the  violent  scatter  occurred  that  was  found  using  both  the  glass 
particles  and  the  HP  295  ball  powder.  As  noted  earlier,  the  only  dif¬ 
ferences  between  the  various  types  of  particles  and  the  experimental  runs 
involving  these  particles  is  the  particle’s  surface  finish. 

5.  3  UNSTEADINESS  IN  THE  BOUNDARY  LAYER  AND  WAKE 

The  interaction  of  the  shock  wave  with  a  particle  is  analogous  to  the 

impulsive  motion  of  a  sphere  in  a  fluid  in  that  the  particle  impulsively 

sees  and  is  acted  upon  by  the  convective  flow  velocity  behind  the  shock 

front.  In  impulsive  motion  a  certain  amount  of  time  is  needed  for  the 

boundary  layer  and  the  wake  to  reach  a  quasi-steady  condition.  In  order 

to  determine  the  time  it  takes  for  the  boundary  layer  to  become  quasi- 

steady,  a  diffusion  time  may  be  defined  as  the  time  required  for  a  sudden 

(19) 

change  to  spread  by  the  process  of  molecular  or  turbulent  diffusion  . 

2 

The  diffusion  time  may  be  taken  as  6  /  v  where  6  is  the  boundary  layer 
thickness  and  v  is  the  kinematic  viscosity.  The  experimental  displace¬ 
ment  thickness  at  the  90°  point  on  a  circular  cylinder  is^^ 


R 

oc 

V 


=  .  8 


If  we  assume  the  following  values  for  R,  U  ,  and  v 
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R  =  3.  28  x  10~4  ft 
1/  =  1.  611  x  10  4  ft*V  sec 
u  =100  ft/ sec 

OC 

6*g()0  =  1,841  x  10"5ft 

The  diffusion  time  may  be  approximated  using  the  displacement  thickness 
for  a  cylinder  and  is 

. 2  fi 
—  =  2. 12  x  10  sec 
v 

The  minimum  time  from  shock  front  passage  over  the  particles  until  the 
first  picture  was  taken  was  always  greater  than  100  microseconds.  The 
diffusion  time  as  computed  above  is  only  two  per  cent  of  this  time  and 
thus  by  the  time  the  experimental  data  is  taken,  the  boundary  layer  should 
be  quasi-steady. 

Experimental  work  on  impulsive  flow  about  cylinders  has  been  done 
by  Schwabe^^  and  very  recently  by  Sarpkaya^^.  Sarpkaya's  experi¬ 
ments  involved  the  impulsive  flow  of  water  about  circular  cylinders.  He 
found  that  the  drag  coefficient  initially  rises  above  the  steady  state  case 
due  to  the  formation  of  the  vortices  and  then  decays  back  to  the  steady 


Then  5 


90° 


becomes 
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state  case  once  the  vortex  flow  has  become  established.  He  found  that 
his  data  correlated  with  a  non-dimensional  time  r.  The  drag  co¬ 
efficient  had  returned  to  the  steady  state  result  when  r  »  28.  r  is  defined 
as 

U  t 

oo  w 

7  ~  R 

where  =  flow  field  velocity 

t  -  time 
w 

R  -  radius  of  cylinder 

For  U  =  100,  R  =  3.  28  x  10  ^  ft,  and  r  =  28  then  t  becomes 

oo  w 

t  =  91. 8  x  10  **  sec 
w 

It  might  be  noted  that  since  the  wake  formation  time  is  an  order  of  magni¬ 
tude  larger  than  the  boundary  layer  formation  time,  then  in  unsteady  flow 
problems,  the  wake  formation  time  dominates.  Since  t  is  less  than  100 
microseconds,  which  is  the  time  between  when  the  shock  passes  the  par¬ 
ticles  to  when  the  first  picture  is  taken,  the  wake  flow  field  has  reached 
a  quasi -steady  condition  before  the  data  is  taken.  Thus  both  the  boundary 
layer  and  the  wake  have  become  quasi-steady  before  any  experimental 
data  is  taken. 
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5.  4  TURBULENCE  IN  THE  FREE  STREAM 


Free  stream  turbulence  can  cause  a  considerable  shift  in  for  a 

(12) 

particular  Re  as  can  be  seen  from  Torobin’s  and  Gauvin's  '  results 
which  appear  in  Fig.  2.  A  literature  survey  yielded  no  information  on 
turbulence  behind  a  shock  wave.  Since  the  shock  waves  are  weak  and  the 
convective  flow  Mach  numbers  are  low,  any  turbulence  which  is  generated 
in  the  boundary  layer  will  not  be  transmitted  into  the  convective  flow. 
Turbulence  will  not  be  generated  by  the  shock  front  since  Schlieren 
photographs  indicated  the  shock  front  was  plane  and  perpendicular  to  the 
walls  of  the  shock  tube.  Further,  the  fact  that  smooth  sapphire  balls 
yielded  data  close  to  the  steady  state  curve  tends  to  substantiate  a 
turbulent>free  flow  field  assumption. 

5.  5  PARTICLE  ROTATION 

Particle  rotation  could  cause  a  shift  in  the  drag  coefficient  of  a 

particle  in  that  on  one  side  of  the  particle  the  separation  point  will  move 

rearward  whereas  on  the  other  side  it  will  move  forward.  This  will 

cause  the  skin  friction  and  the  form  drag  contributions  to  the  drag  co- 

(22) 

efficient  to  change.  Maccolr  '  conducted  some  experiments  to  deter¬ 
mine  the  effects  of  sphere  rotation  around  an  axis  perpendicular  to  the 

4 

flow.  These  experiments  ranged  from  Reynolds  numbers  of  6.  15  x  10 
4 

to  10.  7  x  10  .  Cp  fell  from  .  52  to  .  48  as  the  ratio  of  equatori  peed, 
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U  ,  to  the  relative  velocity,  U  ,  increased  from  0  to  1.  Davies^ ^  jias 

S  K 

conducted  experiments  in  the  same  Reynolds  number  regime  as  Maccoll 

(24) 

and  his  results  agree  with  Maccoll's.  Pasternak  made  observations 

of  freely  moving  spheres  and  suggests  that  the  ratio  of  U  /U  is  on  the 

s 

(25) 

order  of  five  per  cent.  Luthanderv  has  measured  the  drag  coefficient 

of  spheres  rotating  around  an  axis  parallel  to  the  flow  direction  in  the 

region  of  the  critical  Reynolds  number.  Below  the  critical  Reynolds 

number,  a  U  /U  ratio  up  to  two  had  very  little  effect  on  C^. 

s  D 

Garstang^^  and  Drazin^^  studied  analytically  the  effects  of  rotation 
about  the  axis  parallel  to  the  flow  direction  for  the  Stokesian  drag 
regime  and  found  that  the  drag  coefficient  remained  unaffected  by  the 
rotation.  Due  to  the  above  results  that  rotation  does  not  change  the 
drag  coefficient  significantly  and  since  spherical  particles  are  in¬ 
jected  in  the  present  study  in  a  manner  which  does  not  initiate  rota¬ 
tional  motion,  it  is  felt  that  rotational  effects  are  negligible  in  the  present 
study. 

5.6  ACCELERATION  EFFECTS 

A  considerable  amount  of  experimental  work  has  been  done  on  the  ef¬ 
fects  of  acceleration  on  both  sphere  and  cylinder  drag.  Keiirr^,  did 
experiments  with  accelerating  cylinders.  He  found  that  the  drag  coef¬ 
ficient  versus  Reynolds  number  data  could  be  correlated  by  using  an 


75 


acceleration  modulus  defined  as 


He  found  that  if  A  was  on  the  order  of  .  2  or  greater  then  the  effect  of 

(29) 

acceleration  upon  the  drag  coefficient  was  substantial.  Bugliarelkr  ' 

determined  experimentally  the  drag  coefficient  for  accelerating  spheres 

He  found  his  data  also  correlated  with  the  acceleration  modulus,  A  . 

c 

(13) 

Crowe  showed  analytically  that 


where  Cf  =  skin  friction  coefficient.  The  acceleration  moduli  for  the 

-4  -3  / 

present  study  were  of  the  order  of  10  to  10  .  AC^/ becomes 


-  . 35  x  10 


Thus  acceleration  effects  do  not  seem  to  be  important  in  the  present  study 
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VI.  CONCLUSIONS 


Observations  on  the  drag  coefficient  of  small  spherical  particles  in  an 
incompressible,  laminar,  non-reacting,  continuum  flow  regime  are  sum¬ 
marized  below: 

1.  The  drag  coefficient  of  small  spherical  particles  in  laminar,  non¬ 
reacting,  incompressible  flow  regime  was  found  to  be  consistently 
higher  for  a  particular  Reynolds  number  than  the  generally  accepted 
steady  state  value. 

2.  The  CD  of  the  HP  295  ball  powder,  for  MR  <  .  125,  increased 
as  much  as  85  per  cent  and  the  CD  of  glass  beads  increased 
significantly,  for  .  15  <  MR  <  .  30,  over  the  steady  state  value 
respectively. 

3.  The  CD  and  the  scatter  in  the  CD  data  of  the  HP  295  ball  powder 
and  the  glass  beads  substantially  increased  as  the  relative  Mach 
number  increased,  even  though  the  relative  Mach  number  was  still 
in  a  region  which  is  normally  considered  incompressible. 

4.  The  of  smooth  sapphire  balls  fell  closer  to  the  steady 
state  curve  and  did  not  depend  on  the  relative  Mach  number 
nor  did  the  scatter  in  the  CD  data  increase  as  the  relative  Mach 
number  increased. 
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5.  Photomicrographs  of  HP  295  ball  powder,  glass  beads,  and 
sapphire  balls  indicated  that  the  sapphire  balls  were  relatively 
smooth,  the  glass  beads  somewhat  rougher,  and  the  HP  295 
ball  powder  quite  rough. 

It  is  concluded  that  surface  roughness  can  cause  considerable  shift 
in  0^  for  small  spherical  particles  in  an  incompressible,  laminar,  non¬ 
reacting,  continuum  flow  regime.  Since  one  group  of  particles  in  solid 
propellant  rocket  exhaust  is  relatively  rough,  i.e.  ,  those  metalic  particles 
which  are  originally  cast  into  the  solid  propellant,  their  may  be  several 
hundred  per  cent  higher  than  the  steady  state  results.  The  other  particles 
in  the  rocket  exhausts  are  those  which  are  condensed  from  the  com¬ 
bustion  products.  The  of  the  condensed  particles  could  also  deviate 
appreciably  from  the  steady  state  value  depending  upon  the  relative  sur¬ 
face  roughness  and  deformation  of  the  particles.  Thus,  the  use  of  the 
steady  state  curve  for  velocity  lag  calculations  of  particles  in  rocket 
nozzles  possibly  causes  considerable  error. 

This  study  has  only  touched  upon  the  regimes  a  particle  encounters 
in  a  solid  propellant  rocket  exhaust.  Much  work  is  left  to  be  done  on  the 
influence  of  compressibility,  burning,  and  electric  charges  in  the  slip- 
flow,  transition,  and  free  molecular  regimes. 
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APPENDIX 
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DATA  SUMMARY  -  GLASS  BEADS 
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c— 

LO 

o 

00 

tr- 

05 

00 

05 

c- 

05 

c- 

tr- 

CD 

IT- 

CO 

IT- 

c— 

05 

05 

CM 

CD 

LO 

00 

CO 

LO 

LO 

CM 

CM 

CM 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

eg 

eg 

CM 

eg 

CM 

CM 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

05 

eg 

CM 

CO 

c- 

eg 

LO 

LO 

t-H 

05 

CD 

rH 

tH 

o 

r- 

CD 

05 

o 

c- 

CD 

O 

CD 

CM 

LO 

CO 

oo 

IT- 

05 

CD 

05 

C- 

05 

t— * 

o 

05 

tH 

CM 

00 

LO 

CO 

oo 

05 

t- 

r- 

C— 

CM 

LO 

ir— 

00 

CM 

00 

00 

oo 

t— 

oo 

00 

C- 

oo 

tr- 

00 

00 

00 

oo 

05 

00 

00 

oo 

00 

tr- 

t- 

CD 

CD 

CD 

tr— 

CD 

CD 

CD 

t-H 

I 

t-H 

1 

rH 

1 

tH 

t-H 

i 

tH 

1 

tH 

tH 

rH 

rH 

CM 

tH 

tH 

tH 

tH 

tH 

rH 

tH 

tH 

tH 

tH 

rH 

CM 

tH 

t-H 

rH 

CM 

CO 

W 

W 

W 

w 

W 

w 

w 

w 

w 

W 

w 

W 

w 

W 

w 

w 

w 

w 

w 

w 

w 

w 

w 

w 

1 

O 

1 

U 

l 

a 

1 

O 

LO 

c- 

LO 

05 

o 

tH 

CM 

co 

CO 

LO 

00 

o 

CM 

tr— 

00 

o 

tH 

CO 

TT 

05 

CD 

CD 

CD 

CM 

eg 

CM 

co 

CO 

LO 

LO 

lO 

IO 

LO 

LO 

lO 

CD 

CD 

ir- 

tr- 

IT— 

oo 

00 

05 

05 

05 

05 
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TABLE  2  (cont) 


r 


CM 

<D 


PS 


00 

CO 

LO 

02 

CO 

o 

CO 

00 

00 

CO 

o 

t-H 

LO 

CM 

t2~ 

CO 

LO 

02 

o 

02 

LO 

00 

rH 

00 

LO 

LO* 

t-’ 

M<* 

CO* 

CO 

CO* 

CO* 

CM* 

CM* 

02* 

CO* 

rH 

o 

co’ 

CM* 

oo 

1— < 

CM* 

o 

LO 

00 

02* 

CO* 

c— * 

o’ 

CO 

CO* 

CM 

CO 

rH 

co 

CM 

LO 

O 

rH 

C- 

o 

c- 

LO 

o 

t-H 

i-H 

CM 

CO 

LO 

CO 

r-H 

rH 

02 

t— 

o 

02 

02 

lO 

LO 

LO 

LO 

LO 

LO 

CO 

CO 

LO 

LO 

CO 

LO 

LO 

CM 

M1 

Tf 

tJH 

CO 

CO 

LO 

CD  CO  ^  ID  -h 

ID  CO  rH  H  CO 

CD  CD  Is-  h  T* 


10  CO  00  05 

CO  H  CO  ^  CO 

0  CD  CD  CO  CD 


^  O  O  W  03  O 
CD  H  I-  ^  CO  ID 
CD  IT-  <£>  <D  t—  C- 


C—  CM  00  C—  CM 

O  CO  CO  (M  ID 

t2*  C—  £—  c—  c— 


02  Tf  CO  CD 

00  02  O  r~H 

CO  CD  h  h 


CM  lO  CM 
H  W  CO 
h  h  CO 


LO  CO  lO  CO  CO  ID  CO  03  CO  00  CO  02  O  CM  LO  LO  CD  CD  CM  O  h  CO  lO  ^  02 


03  03  O  CD  LD  COl-C^^CM 
^  CD  CO  ID  h 

H  i— (  t-H  rH  H  t — I  rH  t-H  t-H  y—i 


W  CD  CO  H  ID 
h  CD  D  03  O  O 


02  00  02  02  02 

o  o  o  o  o 


CD  LD  O 

o  o  o 


CM  CM 


CM  ^  CM  O  O  CO 

CD  *  •  •  •  • 

q,  ID  03  O  ID  CM 

o  r/3  |>  ID  ID  CM  CD 

\  00  O  02  LD  H 

+->  £2*  CO  C—  t—  00 


C—  OO  00  H  O 

C-*  y—4  CM 

O  ID  CD  03  CM 

CD  02  00  CM  CM 

CO  00  00  c-  02 


O  O  O  O  cm  o 

00  O  CM*  02*  02*  O 
LO  lO  M"1  02  CO  CO 
02  CD  CM  00  ^  O 

CO  D  D  ^  ^  ID 


^  CD  C\]  h  CM  CM  00  CM  00 


O  CD  CO 


LO  LO  OO  CM  M"1 

CM  02  02  O  CO 

02  CM  02  O  CO 

^  ID  ^  ID  LD 


L-  CD  CO  CD 
^  H  CO  ID 
CD  CO  CM  ^ 
^  ^  Tf 


LO  CO  02 
^  t2-  CO 
CO  02 
LO  CO  LO 


X5 


t-H 

S 


0) 

0c2 


Q 

U 


Tf 

o 

o 
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o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

1 

o 

1 

o 

1 

O 

1 

o 

1 

o 
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o 

1 

o 

t-H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

▼H 

t-H 

^H 

t-H 

t-H 

t-H 

t-H 

rH 

t-H 

t-H 

t-H 

t-H 

t-H 

rH 

t-H 

t-H 

t-H 

rH 

t-H 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

LO 

LO 

00 

LO 

t-H 

o 

t-H 

co 

LO 

CO 

LO 

LO 

02 

LO 

o 

o 

LO 

o 

LO 

LO 

O 

O 

LO 

o 

LO 

o 

LO 

LO 

CO 

02 

LO 

CM 

CO 

o 

CD 

LO 

02 

t- 

CO 

lO 

CO 

CO 

CM 

oo 

CO 

LO 

XT 

CO 

rH 

00 

LO 

LO 

CO* 

CO* 

CO* 

CO* 

CO* 

C— 

CO* 

• 

t- 

CO* 

CO* 

CO* 

CO* 

CO* 

^* 

CO* 

co* 

co* 

CO* 

co* 

CO* 

CD* 

CO* 

CD 

CO* 

CO* 

co’ 

CO* 

CO* 

r- 

C— 

t-H 

LO 

LO 

CD 

CO 

CO 

CM 

12- 

12- 

CO 

CO 

rH 

co 

t-H 

CO 

02 

c— 

CO 

o 

LO 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

00 

co 

02 

02 

02 

02 

02 

02 

oo 

02 

00 

oo 

o 

o 

o 

t-H 

rH 

t-H 

t-H 

rH 

t-H 

t-H 

rH 

t-H 

rH 

t-H 

t-H 

t-H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

rH 

t-H 

rH 

02 

02 

LO 

00 

00 

02 

CM 

CM 

CM 

CO 

lO 

rH 

rH 

LO 

02 

oo 

LO 

LO 

LO 

LO 

LO 

oo 

02 

02 

02 

oo 

oo 

00 

CO 

co 

oo 

00 

O 

O 

02 

o 

o 

o 

O 

LO 

LO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

LO 

LO 

CO 

r- 

[r- 

o 

o 

o 

o 

o 

o 

t-H 

t-H 

o 

t-H 

t-H 

t-H 

rH 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

T-H 

t-H 

t-H 

H 

t-H 

t-H 

t-H 

t-H 

t-H 

t-H 

^H 

rH 

tH 

t-H 

rH 

rH 

t-H 

rH 

rH 

t-H 

t-H 

t-H 

t-H 

t-H 

t-H 

rH 

rH 

t-H 

02 

t- 

LO 

CO 

02 

OO 

o 

02 

rH 

CO 

02 

CM 

OO 

CO 

OO 

t-H 

O 

00 

LO 

CO 

CM 

LO 

c- 

o 

c- 

LO* 

J2-’ 

LO* 

CD* 

02* 

CO* 

LO 

LO 

CM* 

00* 

02* 

CO* 

t-H 

o 

OO 

CM* 

02* 

t-H 

OO* 

t-H 

CD 

CO 

02 

o 

CO* 

CM 

oo 

t-H 

CO 

CM 

LO 

O 

t-H 

D- 

o 

c- 

LO 

o 

rH 

t-H 

CM 

CO 

LO 

00 

rH 

rH 

02 

c— 

o 

02 

02 

LO 

LO 

LO 

LO 

LO 

LO 

CO 

CO 

LO 

LO 

CO 

LO 

lO 

CM 

CO 

oo 

lO 

00 

t-H 

00 

c— 

lO 

02 

CM 

rH 

CM 

c- 

oo 

oo 

CO 

CM 

CM 

o 

LO 

rH 

o 

00 

02 

12- 

02 

o 

LO 

02 

LO 

LO 

02 

t-H 

t-H 

oo 

CO 

rH 

LO 

CQ 

t-H 

12- 

LO 

t-H 

oo 

oo 

lO 

00 

02 

o 

CM 

rH 

CM 

00 

CO 

CO 

c- 

£— 

CO 

CO 

CO 

CO 

o 

CO 

t- 

co 

02 

12* 

r- 

r- 

r> 

CO 

CO 

12- 

t- 

12- 

12- 

CO 

rH 

CM 

1 

tH 

| 

t-H 

1 

CM 

i 

rH 

1 

t-H 

t 

CM 

1 

t-H 

1 

rH 

1 

t-H 

t-H 

CM 

H 

rH 

t-H 

t-H 

CM 

oo 

LO 

rH 

t-H 

rH 

CM 

rH 

CM 

tH 

6 

u 

a 

a 

a 

O 

O 

a 

a 

O 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

1 

o 

1 

e> 

6 

6 

1 

O 

1 

(J 

1 

O 

1 

o 

CM 

CM 

lO 

LO 

CO 

oo 

oo 

CO 

rH 

oo 

12- 

12— 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

12- 

02 

t-H 

rH 

c— 

CM 

rH 

rH 

rH 

H 

t-H 

rH 

CM 

CM 

CM 

OO 

00 

oo 

CO 

LO 

CO 

12- 

L— 

c— 

t2" 

D- 

12- 

c— 

00 

00 

co 

00 

02 

83 


TABLE  3 


DATA  SUMMARY  -  GLASS  BEADS 
.  15  <  MD  <  .  30 

K 


RUN 

CD 

Re 

M1 

mr 

d 

a 

9 

UR 

ft 

ft/sec 

ft/sec 

127-1 

.  728 

787.  9 

1.  179 

.  239 

6.  70 

X 

IO"4 

23793. 7 

287.  5 

123-1 

.  591 

894.  7 

1.  208 

.  270 

6.  65 

X 

10-4 

25814.  3 

328.  3 

135-1 

.  614 

811.  5 

1.  199 

.  258 

6.  20 

X 

IO"4 

26337. 3 

311.  9 

130-1 

.  781 

689.  8 

1.  168 

.  225 

6.  30 

X 

10-4 

23590.  2 

269.  9 

126-1 

.  681 

869.  8 

1.  194 

.  255 

6.  80 

X 

IO"4 

25725.  8 

307.  9 

154-1 

.  921 

560.  4 

1.  167 

.  216 

5.  29 

X 

10-4 

30660.  3 

258.  5 

158-1 

.  572 

707.  6 

1.  177 

.  241 

6.  10 

X 

10-4 

20499.  7 

289.  9 

164-1 

.  701 

564.  2 

1.  155 

.  213 

6.  30 

X 

IO'4 

16345.  1 

254.  9 

165-1 

.  857 

629.  8 

1.  163 

.  223 

6.  35 

X 

lO'4 

23180.  4 

268.  2 

166-1 

.  857 

615.  0 

1.  157 

.  214 

6.  55 

X 

lO'4 

20300.  4 

256.  9 

167-1 

.  639 

778.  9 

1.  197 

.  258 

6.  20 

X 

10'4 

26551.  7 

312.  8 

169-1 

.  656 

751.  4 

1.  181 

.  244 

6.  35 

X 

io-4 

23458.  8 

294.  9 

203-2 

.  653 

1378 

1.  130 

.  180 

9.  58 

X 

io-4 

13173.  6 

212.  8 

204-2 

.  686 

1290 

1.  120 

.  166 

9.  83 

X 

IO'4 

11265.  9 

195.  8 

205-2 

.  652 

1577 

1.  138 

.  188 

1.  04 

X 

IO'3 

13456.  0 

222.  9 

225-1 

.  705 

1269 

1.  126 

.  173 

9.  42 

X 

10-4 

13130.  3 

205.  0 

226-1 

.  584 

1331 

1.  135 

.  185 

9.  12 

X 

10-4 

13132.  2 

220.  3 

231-1 

.  594 

1304 

1.  123 

.  170 

9.  88 

X 

IO-4 

10177.  5 

201.  9 

232-1 

.  736 

1150 

1.  105 

.  150 

1.  01 

X 

IO"3 

9242.  1 

177.  7 

238-1 

.  638 

1239 

1.  121 

.  167 

9.  63 

X 

IO’4 

10803.  2 

199.  0 

242-1 

.  719 

1351 

1.  128 

.  175 

9.  88 

X 

10-4 

13196.  3 

208.  2 

244-1 

.  712 

1214 

1.  110 

.  159 

1.  00 

X 

10-4 

10089.  5 

187.  4 

247-1 

.  655 

1158 

1.  118 

.  163 

9.  27 

X 

IO-4 

10735.  1 

192.  5 

250-1 

.  631 

1296 

1.  134 

.  184 

8.  97 

X 

10-4 

14163.  7 

219.  2 

283-1 

.  788 

803.  1 

1.  150 

.  209 

9.  22 

X 

IO'4 

11821.  9 

248.  9 

284-1 

.  834 

608.  8 

1.  119 

.  170 

8.  92 

X 

IO'4 

8071.  0 

200.  9 

286-1 

.  739 

783.0 

1.  147 

.  206 

9.  17 

X 

10-4 

10736.  5 

244.  9 

290-1 

.  524 

896.  9 

1.  163 

.  230 

9.  17 

X 

10-4 

9822.  3 

27  4.  5 

300-1 

.  766 

571.  5 

1.  142 

.  187 

6.  56 

X 

10-4 

14508.  9 

223.  1 

300-2 

.  838 

566.  6 

1.  142 

.  185 

6.  56 

X 

10-4 

15594.  7 

221.  2 

301-1 

.  585 

866.  2 

1.  207 

.  258 

6.  67 

X 

10-4 

23411.  8 

312.  7 

301-3 

.  571 

853.  7 

1.  207 

.  258 

6.  56 

X 

10-4 

23317.  1 

313.  2 

302-1 

.  889 

656.  8 

1.  161 

.  208 

6.  56 

X 

IO-4 

21737. 3 

249.  4 

304-1 

.  629 

804.  1 

1.  193 

.  246 

6.  56 

X 

10-4 

22704.  7 

296.  8 

304-2 

.  657 

774.  6 

1.  194 

.  244 

6.  37 

X 

10  "4 

24082.  6 

294.  6 
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TABLE  4 


DATA  SUMMARY  -  WINCHESTER  WESTERN  HP  295  BALL  POWDER 


RUN 

CD 

Re 

M1 

™R 

d 

a. 

n 

UR 

ft 

ft/sec 

ft/sec 

389-1 

.  939 

596.  9 

1.  101 

.  121 

6.  56 

x  10”4 

11678.  3 

142.  6 

390-1 

.  969 

603.  1 

1.  100 

.  122 

6.  56 

*10-* 

12299.  6 

144.  1 

390-2 

1.  004 

595.  3 

1.  100 

.  124 

6.  56 

*10-* 

12426.  3 

142.  3 

391-1 

.  927 

614.  5 

1.  101 

.  124 

6.  56 

X  10'* 

12227.  2 

146.  8 

577-1 

.  685 

662.  6 

1.  071 

.  099 

9.  32 

-4 

X  10  „ 

5681.  3 

116.  5 

577-2 

.  758 

716.  7 

1.  071 

.  099 

1.01 

X  1°~Q 

5771.  7 

116.  2 

578-1 

.  784 

692.  4 

1.  068 

.  096 

1.  02 

X  10  „ 

5499.  8 

112.  1 

578-2 

.  691 

752.  7 

1.068 

.  097 

1.  09 

X  104 

4618.  3 

113.  5 

581-1 

.  783 

665.  3 

1.  067 

.  094 

9.  92 

x  10* 

5472.  3 

110.  8 

587-1 

.  690 

683.  4 

1.  068 

.  094 

1.  03 

x  10"® 

4596.  2 

110.  0 

591-1 

.  919 

689.  0 

1.  069 

.  096 

9.  78 

xio;* 

6832.  7 

112.  1 

593-1 

.  774 

668.  8 

1.  068 

.  095 

9.  73 

x  10  r. 

5557.  4 

110.  1 

596-1 

.  655 

744.  3 

1.  068 

.  097 

1.  05 

x  10 

4591.  2 

113.  2 

607-2 

.  677 

734.  4 

1.  070 

.  097 

1.  03 

X  10  ^ 

4831.  5 

112.  0 

610-1 

.  701 

724.  0 

1.  070 

.  098 

1.  00 

x  1C' 

5280.  0 

113.  6 

612-1 

.  777 

688.  2 

1.  069 

.  098 

9.  57 

XIO'* 

6060.  1 

113.  2 

613-1 

.  530 

737.  7 

1.  069 

.  100 

1.  00 

X  10  , 

4160.  1 

116.  1 

614-1 

.  692 

722.  8 

1.  069 

.  098 

1.  00 

x  10 

5218.  2 

113.  8 
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TABLE  5 


DATA  SUMMARY  -  SAPPHIRE  BALLS* 


RUN 

CD 

Re 

M1 

“r 

a 

2 

UR 

ft/  sec 

ft/sec 

701-1 

.  581 

1477 

1.  101 

.  149 

3445.  3 

174.  1 

720-1 

.  530 

1462 

1.  100 

.  146 

3048.  7 

170.  8 

723-1 

.  548 

1479 

1.  100 

.  148 

3229. 0 

172.  9 

725-1 

.  572 

1474 

1.  099 

.  148 

3369. 9 

173.  0 

726-1 

.  556 

1487 

1.  101 

.  149 

3334.  6 

174.  4 

728-1 

.  561 

1467 

1.  099 

.  147 

3294.  2 

173.  0 

728-2 

.  571 

1468 

1.  099 

.  147 

3356.  4 

173.  0 

735-1 

.  624 

855.  9 

1.  099 

.  149 

2181.  2 

175.  3 

740-1 

.  473 

956.  0 

1.  107 

.  160 

1958. 2 

187.  8 

741-1 

.  533 

968.  8 

1.  108 

.  162 

2259.  9 

189.  9 

743-1 

.  521 

955.  2 

1.  107 

.  160 

2160.  1 

188.  0 

746-1 

.546 

949.  4 

1.  108 

.  161 

2266.  0 

189.  0 

751-1 

.  465 

1178 

1.  183 

.  257 

4016.  1 

307.  5 

756-1 

.  521 

933.  3 

1.  104 

.  156 

2048.  6 

183.  0 

756-2 

.  531 

932.  9 

1.  104 

.  156 

2083.  8 

182.  9 

758-1 

.  525 

923.  3 

1.  103 

.  154 

2010. 0 

180.  6 

761-1 

.  534 

920.  8 

1.  102 

.  154 

2035.  0 

184.  0 

767-1 

.  501 

1068 

1.  138 

.  202 

2988. 7 

238.  8 

767-2 

.  525 

1069 

1.  138 

.  202 

3134.  0 

238.  9 

775-1 

.  522 

1036 

1.  135 

.  193 

2887. 4 

228.  5 

776-1 

.  484 

1116 

1.  151 

.  217 

3268.  1 

257.  6 

777-1 

.  520 

1088 

1.  140 

.  205 

3209.  1 

242.  5 

778-1 

.  492 

1106 

1.  147 

.  213 

3213.  3 

252.  2 

779-1 

.  520 

1083 

1.  138 

.  202 

3132.  5 

238.  4 

781-1 

.  549 

1061 

1.  136 

.  200 

3199.  2 

235.  8 

783-1 

.  585 

1038 

1.  135 

.  198 

3313.  1 

233.  9 

820-1 

.  488 

1181 

1.  175 

.  250 

4061.  5 

297.  8 

821-1 

.  527 

1158 

1.  177 

.  250 

4301.  4 

297.  5 

828-1 

.  506 

1224 

1.  196 

.  276 

4895. 5 

330.  5 

833-1 

.  492 

1234 

1.  227 

.  313 

5590.  5 

379.  6 

*  -  3 

Diameter,  d,  is  always  1.  3021  x  10  ft  for  sapphire  balls. 
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